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A theory of radar scattering from a rough planetary 
surface based strictly on the geometrical-optics approx- 
imation has been published by Muhleman [1964]. It 
was shown in that paper that the intensity of the 
scattered radiation back to the radar from an element 
of area on a sphere of radius R is (unit flux illumination) 
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where is the angle of incidence (and reflection for the 
radar backscatter case), and p{0)l2n is a probability 
density. The geometrical-optics approximation as- 
sumes that the surface is made up of flat facets whose 
normal vectors are tilted from the mean-surface normal 
by the random angles with a uniform azimuthal dis- 
tribution. Thus, p(6)I2tt is the probability density of 
these tilt angles. 

The density function was normalized to unity over 
the angles in the referenced paper which yielded the 
correct angular spectrum of the scattered energy but 
not the correct total power return. The proper normal- 
ization is obtained by requiring that the total area of 
the facets when projected to the mean spherical sur- 
face by the insertion of cos sum to the area of the 
sphere. Thus the normalized returned intensity is 
found to be 
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where p(0), so normalized, is now identically the angular 
backscatter function. The total intensity back to the 
radar is then obtained by integrating over the hemi- 
sphere, 
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Now it is of major importance in radar astronomy to 
be able to interpret the total measured power in terms 
of the directivity of the surface (related to the back- 
scatter law), and the reflectivity of the planet's surface 
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materials (related to the electrical parameters of the 
surface). If the planet's pulse response or Doppler 
spectrum is known from measurements, then p(0) is 
known and the separation of directivity g and reflec- 
tivity p may be carried out with (3). The directivity 
is defined as the ratio of the backscatter intensity to 
the mean intensity scattered into all directions. The 
mean intensity is clearly ttR 2 /47t; therefore, from (3), 
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1 . The Planetary Scattering Law 

Muhleman [1964] has derived a theoretical scattering 
law based on experimentally distributed height and 
length variations of the scattering facets which fits 
the observed lunar radar scattering measurements 
over a wide range of wavelengths. This law contains 
a single parameter, a, which varies with wavelength. 
It is our purpose here to determine the applicability 
of this law to 12.5-cm radar measurements of Venus 
and Mercury, to determine the statistics of roughness 
for these planets, and to determine the 12.5-cm direc- 
tivities. The unnormalized probability density so 
obtained is 
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[sin + a cos Of ' 

which yields a directivity of 

g^-l + a-2a 2 ; a < 1.0. 

It can be shown that a is the one-dimensional mean 
slope, i.e., the expectation value of tan computed 
with (5). A physically more meaningful statistic is 
the two-dimensional mean slope, where the mean is 
taken over the spherical planetary surface and is de- 
fined by 
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Using (5), we get 

£(tan 6) = Va 



0.09 



(7) 
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2. Applications to Venus and Mercury 

A typical observed radar spectrum of Venus taken 
with the JPL-Goldstone radar system (operated at 
12.5 cm) is shown in figure 1. Also shown by the 
solid line is the spectrum arising from the theoretical 
scattering law with a = 0.10. The agreement is 
excellent. The deviation from the model at a fre- 
quency of 0.5 is due to the presence of a well-under- 
stood spectral feature [Carpenter, 1965]. This value 
of a yields #=1.08, which means that Venus deviates 
very slightly from an isotropic scatterer at 12.5 cm 
(i.e., smooth sphere). From (7), E(tan 6) = 0.316. 

The observed spectrum of Mercury measured with 
the Goldstone system in 1963 is shown in figure 2. 
The quality of the observations is considerably inferior 
to that of the Venus measurements. Consequently, 
two a parameters were used to bracket the measure- 
ments. Figure 3 shows the 1965 Goldstone measure- 
ments of Mercury, which are somewhat smoother, 
primarily due to the use of less resolution. The model 
with a = 0.2 is in good agreement with the data. 
An a =0.2 yields g=1.16 and £(tan 0)=O.45. Thus 
Mercury appears considerably rougher than Venus 
at 12.5 cm and it can be shown that Mercury's rough- 
ness is very similar to that of the Moon. 
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FIGURE 2. Observed spectrum of Mercury and theoretical models 
with a = 0.15 and 0.20. 




MERCURY- MAY 65 
10 c/s WINDOW 

MODEL - a =0.20 

OBS. SPECTRUM 



250 



FREQUENCY c/s 



FIGURE 1. Observed spectrum of Venus and the theoretical scat- 
tering model with a — 0.1. 



FIGURE 3. Observed spectrum of Mercury and model assuming an 
88 rotational period (1965). 



This interpretation of the Mercury spectra assumes 
a synchronous rotation rate for the planet which would 
cause a spectral width of about 70 c/s. The observed 
decay of the spectral wings supports this assumption. 
However, Pettengill [1965J reports a rotational period 
of 59 ± 5 days (as opposed to 88 days for synchronous 
rotation). A critical discussion of the rotational rate 
of Mercury is outside the scope of this paper, but we 
can say that the more rapid rate is not consistent with 
our observations. 
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